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SUMMARY 


A  unified  treatment  of  the  SAR  and  ISAR  imaging  processes  is  presented,  in 
which  each  emerges  as  a  special  case  of  a  more  general  hybrid  SAR/ISAR 
process.  Several  such  special  cases  are  considered,  demonstrating  the 
simplicity  and  flexibility  of  our  analysis  whose  results  are  also  verified 
directly  by  computer  simulation. 
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Introduction 


Any  ship  on  the  sea  surface  is  subject  to  rotational  motions  (yaw,  pitch 
and  roll)  due  to  the  sea’s  movements.  Furthermore  it  is  frequently  convenient 
to  survey  maritime  scenes  from  a  moving  radar  platform.  Both  these  motions 
will  affect  the  imaging  of  a  ship  target  and  can  be  exploited,  through 
suitable  signal  processing,  to  enhance  the  imagery  obtained.  Thus,  if  a 
stationary  target  is  observed  from  a  moving  platform  it  is  possible  to 
synthesise  coherently  a  large  effective  aperture  and  so  to  enhance  azimuthal 
resolution  (SAR).  Conversely  differential  Doppler  processing  of  the  returns 
to  a  stationary  platform  from  a  rocking  ship  can  provide  more  detailed  across 
range  information  on  the  structure  of  the  target  through  the  technique  of 
inverse  synthetic  aperture  radar  (ISAR).  Target  motions  and,  in  particular, 
rotations  will  degrade  SAR  imagery;  ISAR  images  are  also  adversely  affected 
by  radar  platform  motions.  Unfortunately  the  standard  descriptions  of  SAR  and 
ISAR  are  couched  in  terms  which  make  the  analyses  of  these  effects  complicated 
and  their  physical  basis  and  prac:ical  implications  rather  obscure. 
Consequently  it  is  desirable  to  construct  a  theory  of  Hybrid  SAR/ISAR  imaging 
which  treats  target  and  platform  motions  on  an  equal  footing,  both  to  analyse 
the  degradation  of  conventional  SAR  and  ISAR  images  and  to  identify  ways  in 
which  these  hybrid  images  may  be  processed  optimally.  Such  a  theory,  and  some 
of  its  preliminary  applications,  will  be  described  in  this  report. 

In  broad  outline  there  are  three  concepts  which  underpin  our  analysis. 
These  are  the  description  of  the  imaging  process  in  terms  of  Fourier 
transformation,  the  concise  and  correct  representation  of  rotations  and  the 
analysis  of  image  degradation  in  terms  of  the  manner  in  which  spatial  Fourier 
components  of  the  target  are  selected  by  rotational  motions.  The  first  of 
these  ideas,  which  was  introduced  into  the  discussion  of  radar  imaging  by 
Walker^**  allows  us  to  unify  the  description  of  SAR  and  ISAR  and  so  paves  the 
way  to  a  systematic  analysis  of  the  imagery  obtained.  The  description  of  the 
rotational  motions  will  require  the  introduction  of  a  formalism  which,  while 
it  is  perhaps  a  little  unfamiliar  at  first  sight,  provides  us  with  a  powerful 
tool  for  the  analysis  of  the  effects  of  target  and  platform  motions  in  a 
general,  simple  and  readily  visualised  manner.  Finally,  by  bringing  together 
these  basic  ideas  it  is  possible  to  characterise  the  point  spread  or  imaging 
function  in  terms  of  simple  geometrical  concepts  which  greatly  facilitate  the 
discussion  of  the  imaging  process.  The  relationship  between  our  description 
of  the  SAR  and  ISAR  processes  and  the  standard  analyses  is  discussed  in  an 
appendix.  Similarly  the  details  of  the  rotation  formalism  have  been  discussed 
elsewhere^^^  and  will  be  quoted  here,  supported  only  by  plausibility  arguments 
rather  than  by  a  full  analysis.  Having  set  up  our  formalism  we  will  analyse 
the  SAR,  ISAR  and  hybrid  SAR/ISAR  processes  in  some  detail  and  make  a 
preliminary  comparison  with  computer  simulation  results,  which  provide  a  most 
satisfactory  confirmation  of  our  analysis. 
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Imaging  and  Fourier  Transformation 


In  this  section  we  will  discuss  the  way  in  which  the  Fourier  transform 
arises  naturally  in  the  description  of  the  imaging  process  and  can  be  used  to 
characterise  the  images  obtained  in  a  simple  fashion.  This  will  lead  us  to 
a  general  'master  formula'  which,  when  supplemented  by  our  subsequent 
discussion  of  rotational  motions,  will  allow  us  to  analyse  shifts,  distortions 
and  defocussing  effects  in  a  systematic  fashion. 

We  consider  a  target  whose  position  and  orientation  vary  with  time  t  and 
whose  spatial  extent  over  coordinates  ^  is  specified  by  the  scattering 
function  V(r^,t).  The  following  reasonable  assumptions  can  be  made:  That 
there  is  no  shadowing  or,  more  realistically,  that  shadowing  is  unchanged  over 
angles  of  view  sampled  in  the  imaging  process  and  that  there  are  no  non-linear 
or  dispersive  effects.  For  a  monostatic  radar  platform,  positioned  at  r(t) 
(we  include  all  centre  of  mass  relative  motions  of  target  and  platform  in 
r(t)),  illuminating  the  target  with  EM  radiation  of  wavenumber  k,  the  received 
signal  is  proportional  to 


S(k,t)  -  I  V(rQ,t)exp(2ik|rQ  -  r(t) I )d^ro  (1) 

If  the  range  r(t)  is  large  compared  with  the  spatial  extent  of  the  target  then 
this  quantity  is  well  approximated  by 


S(k,t)  -  j  .rQ)d3ro 


(2) 


where  the  scattering  vector  k(t)  is  given  by 


k(t) 


(3) 


and  the  tilde  denotes  Fourier  transformation: 


f  (k) 


-ik.r„  , 

-  -°f(r^)d3ro 


(4) 
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Thus  the  received  signal  can  be  represented  as  the  product  of  a  range 
dependent  phase  factor  and  a  spatial  Fourier  component  of  the  target 
scattering  function.  If  the  range  r(t)  is  known  then  the  phase  factor  may  be 
removed  (this  is  the  essence  of  standard  motion  compensation  procedures)  and 
the  Fourier  component  V(2k(t),t)  determined.  Consequently  successive  returns 
measured  at  several  r(t)  will  sample  the  spatial  Fourier  components  of  the 
target  and  so  provide  us  with  the  raw  material  for  the  imaging  of  the  target 
by  inverse  Fourier  transformation. 

The  precise  way  in  which  target  and  platform  motions  select  these 
Fourier  components  will  be  discussed  in  following  sections  of  the  paper. 
Nonetheless  it  is  worthwhile  at  this  point  to  consider  the  imaging  process  a 
little  further. 

Let  us  suppose  that  we  have  sampled  Fourier  components  f(k)  of  a 
function  f(r)  over  a  limited  set  Z  of  k  vectors.  We  now  attempt  to 
reconstruct  f(r)  (as  f'(r))  by  inverse  Fourier  transformation. 

f'  (r)  -  — r  d^ke^-- -f  (k)H3(k)  (5) 
(2n)3  J 


where 


H3  (k)  -  1  k  G  Z 

-  0  k  €  Z 


Substitution  of  (4)  then  gives  us 


f'  (r)  -  j  H3{r  -  r'  )f (r' )d3r' 


(6) 


where  H3(r  -  r ’ )  is  the  point  spread  or  imaging  function  given  by 


H3{r-r')  -  — r  d^ke^-- -  ^H3(k)  (7) 
(2n)3  J 

which  relates  the  form  of  the  image  we  obtain  to  our  ability  to  sample  k 
space.  Thus,  for  example,  if  we  knew  f(k)  only  in  the  plane  containing 
vectors  k^^  then  we  would  set 


k 


k,  e 


-1  Hj. 
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where  e^^  is  the  unit  vector  perpendicular  to  this  plane  and  set 

^  -  0 

-  0  otherwise 

Consequently 

~  )  "2^  ~  -h'  ^  (S) 

(a  2  dimensional  delta  function)  and 


^  i  f(lH  (9) 

—  cc 

obtaining  an  image  which  is  the  projection  of  f  (r)  onto  the  plane 
containing  the  known  k  vectors. 


This  result  is  of  course  the  well  known  'Fourier  slice'  Theorem,  Should 
our  knowledge  of  f(kH^  restricted  further  in  the  k^  plane 


H3(1Sh  ^  “  H2(kH)  -  0 

-  0  otherwise 


(10) 


then  our  projected  image  will  also  be  defocussed: 

-h’  ■  i  1  «2(rH-rH')  I  dr^'f(rH'  -e^r^')  (11) 

where 

H2(r„-r„')  -  '  f  d2k„  ^  ,  ,, 

For  example  knowledge  of  f(kH)  the  rectangular  domain 
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^01  -  a  <  •^Hl  <  ^01 

^02  -  ^  <  >^H2  <  ^02  * 


would  result  in  a  sine  function  form  for  H2: 


H2(Ih  -  Ih'  ) 


giho-  ^  sina(rHi  -  thi'  )  sinb(rH2  -  '> 

n^Crni  -  thi'  )  (rH2  "  ^H2'  ) 


In  the  general  imaging  problem  k  will  not  be  sampled  in  a  plane  but  on  a  more 
complex  surface  determined  by  the  platform  and  target  motions.  Nonetheless 
it  is  convenient  (and  consequently  is  customary  practise)  to  process  data  (via 
a  F.F.T.  algorithm  or  some  analogue  thereof  (See  Appendix  A))  as  if  it  were. 
This  will  introduce  further  degradation  of  the  imagery,  which  can  be  analysed 
as  follows.  Let  us  take  the  plane  containing  k^  be  tangent  to  the  surface 
containing  the  sampled  k  vectors  such  that  both  surfaces  contain  the  vector 
k^  and  that  centred  around  k^  (c.f  (13)).  (We  will  see  in 
subsequent  sections  that  the  sampled  k  vectors  lie  on  a  cone  like  surface 
generated  by  the  rotation  of  ko  through  relatively  small  angles.  This  readily 
visualised  situation  may  help  motivate  the  following  discussion) .  In 
processing  the  data  as  if  the  k  vectors  were  coplanar  we  perform  a  modified 
Fourier  inversion*  of  the  type 


1 

(2n)3 


t  (}Sh 


(14) 


where  6k (k^)  is  the  vector  correcting  the  approximation  to  the  sampled  k 
vector  to  its  correct  value.  As  the  support  of  H2(kj^)  is  small  and  centred 
on  k^  we  may  expand  6k  in  a  Taylor  series  around  k^; 


6ka  -  1  CqPy  AkHp  AkHy  +  O(Ak^) 


Here 


Ca3Y 


SAkHpSAkHy 


Akn  -  0 


(16) 
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and  we  have  used  the  repeated  suffix  tensor  notation  Thus  we  may  write 


■  -2^  j  ^  »2'  ) 


where 


H2'  (£„,£'  ) 


(2n) 


^ -  J  d^AkH. 


eJ'Pl-'i  ra'CapyAkHpAkny  H2(AkH) 


(17) 


This  very  general  formula  allows  us  to  characterise  the  defocussing  of  the 
image  due  to  quadratic  errors  in  the  approximation  of  sampled  k  vectors  by  a 
planar  set.  In  Appendix  B  we  evaluate  H2'  )  assuming  a  Gaussian  model 

for  H2(Ak^)  [This  is  analytically  tractable  and  yields  readily  interpretable 
results;  a  rectangular  H2  (cf  (13))  yields  results  expressed  in  terms  of 
Fresnel  Integrals  which  are  not  as  readily  visualized]  and  demonstrate  how 
polar  format  processing  effects  arise  as  a  simple  special  case  of  this 
result.  In  the  context  of  SAR  and  ISAR  imaging  we  expect  6k  to  have  a 
component  perpendicular  to  the  k^^  plane,  specified  by  the  curvature  tensor  of 
the  surface  containing  the  sampled  k  vectors.  In  the  following  sections  we 
will  show  how  this  can  be  calculated  using  relatively  simple  geometrical 
ideas.  Furthermore  it  will  be  shown  that  the  corruption  of  the  image  by  errors 
in  the  estimation  of  centre  of  mass  motion  r^  (incorrect  motion  compensation) 
may  be  analysed  quite  straightforwardly  within  the  framework  we  have  just 
presented. 

The  Selection  of  Spatial  Fourier  Components  bv  Target  and  Platform  Motions 

We  now  return  to  Equation  (2)  and  consider  how  target  and  platform 
motions  select  the  spatial  Fourier  components  of  the  target  scattering 
function  VQ(rQ,t)  measured  in  a  SAR  or  ISAR  observation.  We  recall  that  all 
c.m  relative  motion  is  incorporated  in  r(t)  so  that  the  temporal  variation  of 
entirely  to  rotational  motions  of  the  target.  Furthermore, 
as  k(t)  is  of  constant  length  (see  (3)),  its  temporal  variation  can  also  be 
expressed  as  a  rotation.  Consequently,  before  addressing  the  main  topic  of 
this  section,  we  should  discuss  the  representation  of  rotations  we  will  use 
in  our  analysis.  Rotations  of  rigid  bodies  are  commonly  parameterised  in 
terms  of  Euler  s  angles  ,  particularly  in  the  contexts  of  classical 
dynamics^  and  the  quantum  theory  of  angular  momentum  However, 

rotations  can  be  described  in  many  other  ways  (eg.  in  terms  of  Cayley-Klein 
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parameters  and  quaternions^^ ^ )  and  in  this  work  we  have  found  the  axis-angle 
parameterisation  to  be  particularly  useful.  In  this  case  the  rotation  of  a 
rigid  body  in  three  dimensions  is  defined  by  a  unit  vector  (axis)  n,  which 
remains  unchanged  by  the  rotation  and  an  angle  0  which  defines  the  magnitude 
and  direction  of  the  rotation.  This  parameterisation  has  been  discussed 
briefly  in  several  discussions  of  the  quantum  theory  of  angular  momentum 
but  it  is  only  recently  that  its  utility  in  practical  calculations  has  been 
fully  recognised.  A  discussion  of  its  properties  is  given  in  ,  where 
derivations  of  the  results  we  will  use  are  given  and  to  which  the  interested 
reader  is  referred  for  further  details. 

Let  us  consider  a  vector  r,  which  is  rotated  through  an  angle  6  around 
an  axis  n.  r  can  be  resolved  into  components  parallel  and  perpendicular  to 
n,  the  parallel  component  n(n,r)  remaining  unchanged  by  the  rotation.  The 
perpendicular  component  (n  A  r)  A  n  is  transformed  by  the  rotation  into 
(n  A  r)  sin  6  +  (n  A  r)  A  n  cos  0.  Consequently  the  rotation  specified  by  n, 
0  takes  r  into  r'  given  by 


r’  =  n(n  .  r)  +  (n  A  r)  sin  0  +  (n  A  r)  A  n  cos  0 


=  (1  -  cos  0)  n  (n.r)  +  cos  6  r  +  sin  0  n  A  r 


(18) 


which  we  may  express  in  tensor  notation  as 


r  =  M  (n,  0)  .  r 


(19) 


The  Cartesian  components  of  the  dyadic  M(n,  0)  are  readily  identified  as 


-  (1  -  cos0)nQnp  +  cos06a(3  +  sin0  SQyp  ny 


(20) 


where  is  the  Kroenecker  delta,  n^  is  a  component  of  n  and  c^yp  is  the  Levi 
Cevita  or  antisymmetric  tensor  density .  Successive  rotations  characterised 
t>y  (n*! ;  ®l)/(n2/®2)  vector  r  then  give  us  a  vector  r' ’ 


r ' *  =  M  (02/  02^  • 

=  M  (02,02)  •  H  (01 )  .  r  (21  ) 

As  is  demonstrated  in  this  resultant  of  two  rotations  is  equivalent  to 

a  single  rotation  characterised  by  (n,0)  where 
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—  "  ^[(®1®2‘^^(1-ci)(1+C2))ni+siS2n+(1+ci){1-C2) 


(ns2(1-ci)  +  Si  ( 1  -  C2 ))  n2  A 


-  (S1S2+Q  (1  -c)(1  +C2))^  -f  (siS2n+  (1  +ci)(1  -C2))" 


2!1(siS2  +  C^(1  +C2)(1  -  c)|(siS2n  +  (1  +  ci)(1  -C2)) 


+  (1  -  n^)(ns2  (1  -Cl )  +  si  (1  -C2))^ 

In  these  formulae  =  cos  0;  ,  =  sin  6^  (i  =  1,2)  and  Q  =  n).n2,  the 
cosine  of  the  angle  between  the  axes  of  the  successive  rotations.  The  most 
striking  feature  of  these  results  is  that  the  resultant  of  successive 
rotations  about  non  parallel  axes  (n^  n2 )  depends  on  the  order  in  which  the 
rotations  are  performed,  that  is  that  rotations  in  general  do  not  commute 

M  (n2,  02)  .  M  (n^ ,  0^)  #  M  (n^ ,  0^)  .  M  (n2,  02)  (25) 

It  is  also  possible  to  derive  an  equation  of  motion  for  n,  0  describing  the 
rotation  of  a  body  subject  to  a  time  varying  angular  velocity  ^(t).  As  is 
discussed  in  it  can  be  shown  that  the  vector  g(t)  =  0{t)  n(t)  satisfies 
the  differential  equation 


where  is  the  unit  dyadic.  For  small  q  (small  rotations)  we  have  the 
limiting  form 
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^  -  CO  +  ^(wA  3)  (27) 

at  2 

and  note  that  this  differential  equation  is  well  behaved  in  the  vicinity 
of  g  =  0,  which  represents  a  zero  nett  rotation. 

Having  reviewed  our  rotation  formalism  we  turn  to  our  analysis 
of  V(2k(t),t)  .  As  V(rQ;t)  changes  only  as  a  result  of  rotations  of  the 
target  we  may  write  this  in  terms  of  Vq(^)  which  specifies  the  spatial  extent 
of  the  ship  in  the  orientation  adopted  at  time  t  =  0,  as 


V(^,t)  =  Vq  (M  (ngCt),  eo(t))’^  . 


Vo 


where  ng(t),  0s(t)  specify  the  rotation  performed  by  the  ship  in  time  t.  Thus 
we  have 


r  /  -1  \  -2ik(t )  .r* 

V(2k(t),t)  -  I  d^roVo|M^(t)  '  ^r^je  ^ 


I  d3r. 


Vo(ro)  e 


-2ik(t).Mg(t).rQ 


I  Vo(rQ)exp  -2i|Mg(t)  ^•k(t)|.rQ' 


Here  we  have  made  the  change  in  variables 


a  transformation  which  has  unit  Jacobian,  and  have  exploited  the  orthogonality 
property  of  M 


M  -  M‘ 


where  M,M  ^  denote  the  transpose  and  inverse  of  M  respectively.  Finally  we 
note  that,  in  accordance  with  our  earlier  remarks,  we  may  set 
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k(t) 


”p(t).ko 


k 

F(^ 


SO  that 


V(2k(t),t)  -  Vo(2Mg(t)-'' .Mp(t) 


(30) 


showing  explicitly  how  the  target  and  platform  motions  select  out  spatial 
fourier  components  of  the  target  scattering  function  in  a  fixed  orientation. 

ISAR  and  SAR/ISAR  Imaging 

Having  shown  how  the  target  and  platform  motions  select  spatial 
frequency  components  of  the  target  scattering  function  we  will  now  apply  our 
’master  formula’  (17)  to  characterise  SAR,  ISAR  and  hybrid  SAR/lSAR  images  in 
terms  of  these  motions.  Let  us  consider  firstly  the  case  where 


M^(t)  -  M^(t)  ^Mp(t) 


represents  a  rotation  around  a  fixed  axis  n.  The  special  cases  of  SAR  and 
ISAR  with  uniform  platform  and  target  rotations  will  be  covered  by  this 
analysis.  Subsequently  we  will  analyse  the  effect  of  changes  in  the  direction 
of  the  axis  of  rotation  during  the  imaging  process  and  will  then  be  able  to 
address  the  general  SAR  and  ISAR  problems,  as  well  as  the  hybrid  SAR/lSAR 
problem.  We  will  also  consider  the  effects  on  the  imagery  of  errors  in  the 
estimation  of  the  centre  of  mass  relative  motion  r(t)  and  so  be  able  to 
quantify  the  effects  of  incorrect  motion  compensation.  In  this  way  we  hope 
to  show  that  all  the  essential  features  of  imaging  exploiting  target  and 
platform  motions  can  be  incorporated  in  a  relatively  simple  and  systematic 
conceptual  and  computational  framework. 

When  a  vector  is  rotated  about  an  axis  n  the  resulting  vector  k  lies 
on  a  cone  whose  axis  is  n.  This  basic  geometrical  construction  is  shown  in 
Figure  1  .  This  allows  us  to  identify  the  surface  on  which  the  sampled  k 
vectors  lie.  When  we  now  process  data  obtained  in  this  way  by  Fourier 
inversion  (cf(14))  we  will  obtain  an  image  projected  onto  the  plane  tangent 
to  this  cone  and  containing  k^.  Simple  vector  geometry  shows  that  this  plane 
is  normal  to  the  vector  k^  A  (nAk^)  .  The  pertinent  features  of  the 
geometry  of  a  cone  lying  on  a  tangent  plane  are  discussed  in  Appendix  C  while 
in  Appendix  B  we  evaluate  the  tensor  which  we  then  substitute  into  our 
master  formula  (17).  is  now  resolved  into  components  parallel  (k». )  and 
perpendicular  (kj^2)  to  k^  and  H2(kH)  is  assumed  to  have  the  Gaussian  form 
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H2(kH) 


expj  - 

2n  det  D”^]^ 


where  D  is  diagonal 


E  -  (; 


’di  0 

0  d2 


As  is  demonstrated  in  Appendix  B  this  form  of  D  leads  to  an  intensity  point 
spread  function  with  a  covariance  matrix  given  by 


C22  C11 


^  fC22  C11 


^22  ^  <=12 
”3T  di 


where  c^^^  are  given  by 


r2/ko 


l^Q  l^Q 


Q  being  the  angle  between  n  and  k^. 

When  we  identify  n  with  the  axis  of  rotation  of  the  target  and  the 
platform  is  assumed  to  be  stationary  (33)  characterises  the  defocussing  of  an 
ISAR  image.  Conversely  for  a  stationary  target  and  a  platform  whose  motions 
are  generated  by  rotation  around  an  axis  n  (33)  characterises  a  SAR  image. 

Let  us  now  consider  effects  of  the  incorrect  estimation  of  centre  of 
mass  relative  motions.  To  do  this  we  return  to  our  starting  point  (2)  and  set 


Lit)  =  r^it)  +  6r(t) 


where  r^it)  is  our  estimate  of  the  centre  of  mass  motion.  If,  as  is 
reasonable  in  practise,  we  assume  that 

|re(t)  I  >  >  |ro(r) |  ,  |6r(t) | 

we  find  that 


k|£o  ~  £(t)|  =  krg(t)  +  6r(t)  .  k^(t)  -  .  k^Ct) 

where 


ke(t)  -  k 


£e(t) 


^£(t)  may  now  be  decomposed  into  a  gross  error  6r  and  time  varying  parts  due 
to  errors  in  the  estimation  of  centre  of  mass  velocity  and  acceleration 
respectively 


6r(t)  =  6r  +  vt  +  J,  at^ 

Substitution  into  (1)  then  give  us 


where  Mpj(t)  =  MgCt)  ^.Mp(t). 

Consequently  Fourier  inversion  of  data  sampled  with  an  incorrect  estimate  of 
r(t)  will  yield  an  image  of  the  form: 
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r  (rH)  -  ^  J  d3r'  Vo(r'  ) 


iiSo-  -  £h'  - 


I  d2kH( 


+  ikH .  +  6L'  + 


v.kQ(nAko) 
- ^  -  Tu  ) 


0)  I  n  Ak^ 


ci)|  n Ak^ 


X  e  2  r^-  CapykHal^HPH 


^2 


where  we  have  used 


t  =  Kh  •  (n  A  k^)/a)|n  A  k^  I  ^  +  0(k„)^ 


The  effect  of  the  linear  phase  shift  within  the  integral  over  k^  is  to 
shift  the  argument  of  the  imaging  function 


Ih  -  £h'  ^  Ih  ^ 


v.k^in  A  k^) 

wh'TTk^P 


inducing  a  shift  in  the  image  of  the  target,  proportional  to  the  radial 
component  of  v  and  in  the  direction  of  n  A  In  the  standard  SAR  case  this 
shift  is  in  the  'along  track'  direction  while  in  the  ISAR  case  the  familiar 


Cl)  6  r  a  V 
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relationship  is  in  evidence.  The  quadratic  phase  correction  terms  induce  a 
further  defocussing  of  the  image;  unlike  the  'curvature*  defocussing  we  have 
just  discussed  this  defocussing  depends  only  on  the  experimental  geometry  and 
the  errors  in  c.m  estimation  v,  a  and  not  on  the  spatial  extent  of  the  target. 
Reference  to  Appendix  B  allows  us  to  obtain  a  quantitative  estimate  of  these 
defocussing  effects  through  (17). 

We  now  consider  the  effect  of  changes  in  the  direction  of  the  axis 
generating  the  rotation  of  the  k  vectors  within  the  time  during  which  we  form 
the  image.  This  will  bring  about  a  further  contribution  to  6k{kj^)  which  will 
not  in  general  be  perpendicular  to  the  imaging  plane.  In  our  discussion  of 
defocussing  effects  we  need  only  consider  contributions  to  6k  which  are  of 
order  (Akj^)  .  From  (29)  we  see  that  g  =  n0  satisfies  the  differential 
equation 


^  +  1  a)(t)  A  q  (38) 


in  the  limit  of  small  q  (sufficient  for  our  purposes)  where  the  angular 
velocity  varies  in  direction  as  well  as  magnitude.  Thus  we  write 


qit)  =  a)Q(t)  n  +  q^{t) 


(39) 


where 


n  .  0)^  (t )  =  0 

and  n  is  a  fixed  direction. 
Finally  we  set 

t 

eo(t)  -  j  c^(t' )  dt' 

O 

and  write 

g(t)  = 

=  <53{t) 
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Substitution  into  (29)  then  gives 


^  6g(t)  -  w(t)  +  1  ^  ^3(t)  +  (o(t)  A  [go(t)  +  6g(t)j 


which  has,  to  lowest  order  in  u-|(t),  the  solution 
t  t 

6g(t)  -  I  coi(t')dt'  j  ui(t' )  Ag^tt' )dt' 

o  o 


so  that 


t  t 

g(t)  -  eo(t)n  +  I  coi  (t'  )dt'  +  ^  |  eo(t' ) (t'  )dt'  A  n 

We  note  that  the  second  and  third  terms  in  (45)  are  orthogonal 
Consequently,  if  we  retain  only  those  terms  which  are  linear  in  o-j(t) 
that 


e(t)  -  0o  (t). 


n(t)  -  n  + 


1 

^^TtT 


t 

j  wi  (t'  )dt' 

O 


1 


t 

j  0o  (t'  )  co^  ( t'  )  dt'  A  n 
o 


=  n  +  6n ( t ) . 


Thus,  at  some  time  t,  the  vector  will  have  been  rotated  into  k’(t) 
by 


k'  (t)  “  (1  -  cos0o(t))  (n+  6n)  .  ko  ^0  cos0o(t  Ik^ 


+  sin0Q(t)  (n  +  6n)  A  k^ 


The  corresponding  vector  k(t)  obtained  by  a  rotation  of  k^  through  0Q(t 
the  unperturbed  axis  is 


k{t)  “  -  cos0o(t)j  (n.k^)  n  +  cos0o(t)kQ  +  sin0o(t)n  A  k^ 


(40) 


(41) 


(42) 

to  n, 
we  see 


given 


(43) 


)  about 


(44) 
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The  change  in  k  due  to  the  perturbation  of  the  axis  of  rotation  is  given  by 


k'(t)-k(t)  -  [^  -  coseo(t)](  {6n.kQ)n  {n.kQ)5nl 


+  sin0o(t )  5n  A  k 


with  which  we  must  supplement  6k  to  identify  the  corresponding  defocussing 
effects.  We  must  now  evaluate  6r(t)  in  terms  of  G^Ct);  to  do  this  we  must 
postulate  forms  of  w^Ct),  For  short  times  we  set 


w-j  ( t )  =  (i)  ^  t 


so  that 


Ig^w^  a  n 


From  Appendix  C  we  see  that 


A  k^) 
In  A  k^l^ 


so  that  ultimately  we  have 


k'(t)  -k(t)  -  J^L 

^  hTTO 


“1  ^  Ho 


with  which  we  supplement  (17)  to  yield 
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Comparison  with  the  analysis  in  Appendix  B  shows  that  the  effect  of  variation 
in  the  direction  of  the  axis  of  rotation  merely  supplements  the  C22  element 
of  C  by  a  factor 

0)1  E'  •  ( “1  A 
"  4  lnAkj2 

Hybrid  SAR/lSAR  in  which  both  the  target  and  platform  have  motions  provides 
an  interesting  case  in  which  the  direction  of  the  axis  of  rotation  varies  with 
time.  Rotations  about  non-parallel  axes  do  not  commute  and,  as  a  consequence 
of  this,  the  direction  of  the  axis  specifying  the  resultant  of  two  such 
rotations  will  vary  in  time.  Specifically,  if  the  rotation  of  the  target  is 
specified  by  axis,  angle  variables 
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while 


0p(t),  np(t) 


specify  the  rotation  of  the  platform  then  6(t),  n(t)  specifying  the  resultant 
rotation  satisfy 


n(t)  e(t)  -  -6s(t)  n(t)  +  6p(t)np(t)  -  ^  es(t)ep(t)  ^^(t)  Anp(t)  +  OCe^) 


(51 ) 


If  now  we  set 


e5(t)  Sgtt)  =  Qg  cogt 
6p(t)  np(t)  =  ^  u>pt 


then 


e(t)  n  (t) 


%> 


0(t3) 


(52) 


comparison  with  (39)  then  allows  us  to  make  the  identifications 


“o  ■  -  “s“s 

w-|  -  A  ^ 


noting  that  o-]  is  indeed  perpendicular  to  Oq,  as  required  in  our  analysis. 
Consequently  we  see  that  the  non  commutativity  of  rotations  leads  to  an  extra 
source  of  defocussing  in  hybrid  SAR/ISAR  imaging,  which,  as  we  have  seen,  can 
be  accommodated  quite  straightforwardly  in  our  formal  framework. 
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Point  Target  Simulations 


In  this  section  an  initial  verification  of  the  theory  from  the  previous 
sections  is  provided  by  comparing  the  effect  of  mismatched  processing  on  the 
predicted  and  actual  sizes  and  shapes  of  the  image  of  a  point  target.  The 
mismatch  introduced  is  simply  that  which  results  in  normal  practise  from  a 
lack  of  knowledge  of  the  target  rotation  rate  and  from  the  requirement  to 
produce  a  two  dimensional  image  from  three  dimensional  k  vector  samples. 

Two  examples  are  chosen  as  indicated  in  Figure  2,  In  both  cases  a 
target  is  viewed  at  a  low  grazing  angle  by  a  stationary  radar.  The  target  is 
rotated  with  a  constant  rotation  vector  w.  In  case  one  (Figure  2a)  o  is 
vertical,  corresponding  to  a  yaw  motion  of  a  ship.  In  case  two  (Figure  2b) 
0)  is  45®  from  the  vertical,  in  the  plane  containing  the  radar  line  of  sight 
and  the  vertical  direction.  A  number  of  individual  point  scatterers  are 
positioned  on  the  target  in  each  case,  as  shown  in  Figure  2.  The  radar 
parameters  adopted  in  the  simulations  are  listed  below  the  figure. 

It  is  assumed  that  the  radar  has  a  wide  instantaneous  bandwidth,  so  that 
a  range  profile  may  be  derived  for  each  angular  position  of  the  target.  The 
rotation  vector  is  assumed  to  be  unknown,  so  that  range  and  Doppler  migration 
cannot  be  corrected  effectively  prior  to  image  formation.  The  image  is 
therefore  derived  by  simple  Fourier  transformation  of  a  time  sequence  of  radar 
returns  in  each  range  cell  of  the  range  profile.  The  processing  is  (as  is 
shown  in  Appendix  A)  analogous  to  assuming  that  the  measured  k  vector  samples 
come  from  a  rectilinear  two  dimensional  grid  on  a  plane  in  k  space.  As  we 
have  seen  in  our  previous  discussion  k  vectors  are  sampled  from  the  surface 
of  a  cone  (providing  the  rotation  vector  w  does  not  change  direction  during 
the  imaging  process)  and  the  projection  plane  is  tangent  to  this  cone.  In 
Appendix  B  it  is  shown  that  the  errors  6k  in  the  rectangular  grid  processing 
described  above  lead  to  a  distortion  of  the  target  point  spread  function  such 
that 


2 

<x^> 


di 


r2'2 


(53) 


2 

<y^  > 


^2 


(ri '  -  r3'  cotD)^ 

^^2 


r2'" 

Ml 


(54) 


<xy  >  -  - 


r2'  (ri '  “  r3'  cotn) 


(55) 
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where  r-j*,  ^2  ' »  ^3*  position  coordinates  of  the  scatterer  relative  to 
the  rotation  origin  and  d-|,d2/  ^  and  are  as  defined  in  Appendix  B.  Figures 
3,4;  5  and  6  show  the  point  spread  functions  (PSF)  derived  by  simulation  of  the 
rotation  and  radar  processing  of  the  various  scatterers  in  the  two  cases 
illustrated  in  Figure  2.  Table  1  compares  the  moments  of  the  PSF  with  (53), 
(54)  and  (55).  The  excellent  agreement  obtained  indicates  that  the 
approximations  made  in  our  analysis  are  sufficient  to  provide  an  accurate 
representation  of  the  imaging  process. 

It  is  worth  considering  whether  the  results  obtained  from  the  formal 
analysis  used  here  can  be  predicted  by  the  ad  hoc  methods  customarily  employed 
in  the  description  of  ISAR.  Figure  3  shows  the  PSF  for  a  target  at  the  origin 
in  case  1  .  A  scatterer  offset  in  x  (along  the  line  of  sight  (LOS))  undergoes 
Doppler  migration  in  Range  Doppler  processing;  this  explains  the  additional 
spread  in  y  in  figure  4(b).  A  scatterer  offset  in  y  (perpendicular  to  both 
LOS  and  the  rotation  vector)  undergoes  Range  migration,  which  both  produces 
a  spread  in  x  and  reduces  the  time  in  each  range  cell.  This  in  turn  produces 
a  spread  in  the  Doppler  domain  y  and  so  explains  Figure  4(a).  The  spreading 
of  a  point  offset  in  both  x  and  y  is  more  difficult,  if  not  impossible,  to 
derive  (or  rationalise  post  hoc)  in  this  way. 

In  case  2  the  rotation  vector  0  is  not  perpendicular  to  the  LOS.  The 
usual  ISAR  analysis  of  this  case  considers  only  the  projection  of  q  onto  the 
plane  perpendicular  to  LOS.  This  ad  hoc  method  correctly  predicts  the  image 
resolution  of  points  offset  in  directions  perpendicular  to  the  projection  of 
03  (x  or  y  in  our  figure)  but  fails  to  predict  a  mismatch  for  scatterers  offset 
along  the  direction  of  projection  of  o)  (ie.  in  2).  A  mismatch  does  occur, 
however,  as  illustrated  in  Figure  4(d)  and  is  correctly  described  by  the 
formal  methods  used  here.  This  simple  example  demonstrates  the  dangers  of 
analysing  ISAR  from  the  starting  point  of  Range-Doppler  processing,  subject 
to  progressive  corrections  applied  ad  hoc,  and  the  ease  with  which  our 
formalism  encompasses  cases  for  which  the  conventional  analysis  is  either  too 
complex  or  incorrect. 

Conclusions 


In  this  report  we  have  introduced  a  theoretical  framework  within  which 
it  is  possible  to  analyse  the  effects  of  simultaneous  target  and  platform 
motions  on  SAR  and  ISAR  processing  techniques,  by  unifying  both  within  the 
hybrid  SAR/lSAR  process.  Although  our  formalism  is  general  and  of  widespread 
potential  utility,  we  have  shown  that  it  is  readily  applied  to  practical  SAR 
and  ISAR  situations  which,  from  the  standard  Range-Doppler  imaging  standpoint, 
are  difficult  if  not  impossible  to  analyse  correctly.  Two  such  special  cases 
have  been  simulated  on  the  computer,  with  results  which  are  in  excellent 
agreement  with  the  predictions  of  our  analysis. 

While  the  work  reported  here  is  necessarily  of  a  preliminary  nature  we 
feel  that  the  approach  we  have  developed  has  already  demonstrated  its 
practical  applicability  and  that  it  has  the  potential  to  make  tractable  many 
unresolved  problems  occurring  in  the  imaging  of  maritime  targets. 
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APPENDIX  A 


In  this  appendix  we  establish  the  connection  between  standard  SAR  and 
ISAR  processing  and  our  k  space  imaging  formalism.  The  conventional 
description  of  ISAR  processing  starts  with  a  high  range  resolution  radar  of 
bandwidth  B  recording  a  set  of  range  profiles  with  range  resolution  =  2B/c. 
If  the  target  is  rotating  each  scatterer  on  the  target  is  moving  with  respect 
to  the  radar  with  a  range  rate  v^  proportional  to  its  distance  from  the 
target’s  centre  of  rotation,  measured  perpendicular  to  the  radar  line  of  sight 
r^  and  the  rotation  vector  w  ie. 


Vr  -  (coAr).r^  -  A  w)  .r 


(A.1) 


where  r  is  the  position  vector  of  the  scatterer.  The  Doppler  frequency  f^  of 
the  scatterer  is  given  by 


f^  -  -  .^£2(f  Aa)).r  (A.  2) 

Thus,  for  a  small  rotation  angle,  a  frequency  analysis  of  the  signal  in  each 
range  cell  produces  an  output  at  each  frequency  f  corresponding  to  all 
scatterers  at  a  distance  fc/(2fo£o  A  oj)  from  the  centre  of  rotation,  measured 
in  the  direction  (r^  A  w)  .  The  across  range  resolution  is  given  by 


^a 


c 

2Tf  o^o  A  c^t ) 


(A. 3) 


On  converting  this  description  to  spatial  frequency  k  we  find  that,  for  small 
angles  of  rotation,  V(k)  is  measured  in  a  plane  perpendicular 
to  (wAr^)  Ar^  over  an  interval  6kj.  in  the  direction  r^  and  6K^  in  the 
direction  w  A  r ^  where 


6kr  - 


4nB 

c 


6ka 


^oto 


A  ^t 


(A. 4) 


From  the  Fourier  transform  relation  6K  =  2n/Q,  we  find  that  the  equations  for 
^r'  ^a  equivalent  to  those  for  6K^,  6K^. 

For  SAR  range  resolution  is  identical  to  ISAR.  Across  range  resolution 
is,  however,  derived  by  integrating  the  returns  from  a  sideways  looking  radar 
mounted  on  a  moving  platform.  For  an  airborne  radar,  travelling  at  velocity 
V  perpendicular  to  the  radar  line  of  sight  at  a  range  R  from  the  target,  the 
rate  of  change  of  angle  of  view  G  at  the  target  is  given  by 


A1 


d0 

dt 


V 

R 


(A. 5) 


The  usual  method  of  processing  SAR  data  is  to  correlate  the  received  signal 
from  each  range  cell  with  a  linear  FM  waveform.  This  is  equivalent  to  firstly 
correcting  the  phase  of  the  signal  for  the  change  of  range  between  the  radar 
and  secondly  performing  a  Fourier  transform.  By  comparing  the  Fourier 
transform  to  the  frequency  analysis  of  ISAR  processing  and  equating  ^  (A. 5) 
to  w  in  ISAR,  one  shows  that  the  SAR  and  ISAR  processes  are  essentially  the 
same,  the  only  difference  being  that  in  ISAR  the  rotation  is  performed  by  the 
target  while  in  SAR  it  comes  from  the  platform  motion.  The  connection  with 
k  space  processing  for  SAR  is  therefore  quite  analogous  with  that  for  ISAR 
discussed  above. 
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APPENDIX  B 


In  this  appendix  we  evaluate  the  master  formula  (1 
represented  by  the  Gaussian  function 


expj-  ^  (kjj  -  )^)  .D.  (kjj  -  )^)| 

2n[detD"^ 


where 


D  -  (^11  ‘^12] 


[<^21  ^22) 


and  that  (c.f  17) 


6k. r'  -  ^  (kH  -  ko) 


where 


cil  C12 
C21  C22 

Consequently 

H2'  (5;h'£'  ^ 


^iko- ■  £h’  ^ 

(2n)%etD"''|5 


n  iktr.  ( r T_T  —  rTj  ) 

d2kH  e 


1  c 


(D 


i  C)  . kjj 


We  now  focus  our  attention  on  the  integral 


d2kH 


ikn-^ 


®^P  H'2'— H 


B1 


7 )  assuming 


(B.1) 


(B.2) 


(B.3) 


(B.4) 


(B.5) 


where  R  -  r^j  -  rj^'  and  Z  -  D  +  iC 


This  has  the  value 


(detZ)3 


expj-  -^R.Z  ^  .R 


(B.6) 


where 


Z’’'  -  Adj  _ i_ 

=  (detZ) 


In  the  special  case  of  2  x  2  matrices  the  following  identities  hold: 


Adj(D  +  iC)  -  AdjD  +  iAdjC  , 


det{D  +  iC)  -  detD  -  det  C  +  i  D  g  C 


where 


DoC  =  C22dl  1  +  ci  id22  -  ci2d2i  -  C2idi2 


D. Adj  C  +  C. Adj  D  - 


i) 


Consequently,  if  we  set 


detD  -  a,  detC  -  p,  D  g  C  -  y 


(B. 10) 


we  may  write 
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Ad  j  D  +  i  Ad  j  C 

(a  -  3  +  i  Y) 


((a  -  3)  Adj  D  +  Y  Adj  C)  +  i  ((a  -  3)  Adj  C  -  Y  Adj  D) 
((a  -  3)2  +  Y^) 


(B.11) 


SO  that 


j  d^kne^^^'^expj-  ^  k^.Z-k^l 


2n 


[(a  -  3)  +  Y^]" 


exp 


i  tan-1 


Y 

a  -  3 


X  exp 


(a  -  3)  R.Adj  D.R  +  YE-Adj  C.R 


[{a  -  3)^  +  Y^. 


X  exp 


1  (a  -  3)R.AdjC.R  -  YE-AdjD.R 

2  [(a  -  3)^  +  Y^  ] 


(B.12) 


where  we  have  explicitly  separated  the  expression  into  amplitude  and  phase 
factors.  While  the  phase  term  may  be  of  importance  in  future  work, 
particularly  with  reference  to  the  interference  patterns  which  occur  in 
defocussed  images  of  complex  targets,  we  will  restrict  our  attention  here  to 
a  more  detailed  discussion  of  the  amplitude  factor.  To  simplify  the  formulae 
we  make  the  assumption  that  D  is  diagonal  (the  1  or  x  direction  being  along 
k^,  the  2  or  y  direction  being  along  n  A  k^  ) 


ie 


D 


^  di  0  ^ 
.  0  ‘^2. 


B3 


r 


while  C  is  symmetric  and  takes  the  form 


C 


'•=11  ci2  ' 
,ci2  C22^ 


a,  P,  and  y  can  now  be  evaluated  as 


a  =  d^  d2 

P  =  ^22  1  "  *=12 

Y  =  dll  ^22  +  ^22  ^11 


The  amplitude  point  spread  function  now  takes  the  form  of  a  Gaussian 

exp  I  -  ^  R  .  A  .  R  j  (0  13) 

where 


A 


'  (a-(3)  d2  +  YC22 

(a-P)2  + 


-YC12 
(a-P)^  + 


-TCI  2 

(a-p)^  +  Y^ 


(a-P)  d-j  +  YCi  1 

(a-p)2  +  Y^ 


{B.14) 


The  covariance  matrix  of  the  point  spread  function,  which  gives  a  quantitative 
measure  of  the  defocussing,  is  given  by  the  inverse  of  A,  which  takes  the 
relatively  simple  form 
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This  result  forms  the  basis  our  interpretation  of  the  simulation  results  in 
the  body  of  the  text.  It  now  remains  to  evaluate  the  C  matrix  in  a  couple  of 
special  cases  of  practical  importance.  To  make  contact  with  our  simulation 
work  we  consider  the  set  of  k  vectors  generated  by  the  rotation  of  a  vector 
k  around  an  axis  n.  The  magnitude  of  k  varies  with  the  frequency  cycle  of  the 
transmitted  signal.  Imaging  is  carried  out  by  Fourier  transformation  in  the 
plane  normal  to  k^  A  (nA)^)  ,  A  given  k  vector,  of  magnitude  k^  and 
direction  obtained  by  a  rotation  of  y  about  n  from  that  of  k^,  can  be  resolved 
into  components  (1,2,3)  parallel  to  (ko' S  A  k^  A  (nA)^))  as 


k  -  (  kt  (cos^  +  cos  0  sin^  ,  k^  sin^^  sin0,  k^  (  1  -  cos0 )  cos  Q  sin fi) ,  (B.  1  6) 

n  .  k 

where  cosQ  -  _ _ 

^o 

This  vector  is  approximated  in  the  transform  plane  by 


kn  -  (k^,  0koSinQ,  0) 
so  that  we  may  write 

6k  -  (  kt  sin^  Q  (cos0  -  1  )  ,  (k^  sin  0  -  k^  0  )  sinQ, 
kt(1  -COS0)  sinQ  cosQ  ) 


(B.17) 


(B.18) 


Elements  of  the  tensor  may  now  be  obtained  by  differentiation  with 

respect  to  (  k^  ^  1  ”  ^t  ^2  “  6koSinr^  and  setting  k^  =  k^,  0  =  0.  In 

this  way  we  evaluate  the  C  matrix  as 


0 


^2 


^o 


/  t 

-r^  r^cotn 


(B.  19) 


remembering  that  Q  is  the  angle  between  and  n,  so  that  the  covariance 
matrix  takes  the  form 
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The  result  appropriate  to  polar  format  defocussing  is  obtained  by  setting 
cot  n  =  0. 
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APPENDIX  C 


In  this  appendix  we  will  examine  those  features  of  the  geometry  of  a 
cone  lying  on  a  tangent  plane  which  are  required  in  our  analysis  of  SAR  and 
ISAR  imaging.  In  particular  we  will  show  how  to  relate  the  component  of  a 
vector  lying  on  the  cone  resolved  onto  the  tangent  plane  to  its  component 
perpendicular  to  the  plane  and  so  derive  an  expression  for  the  curvature  of 
the  cone.  Sufficient  detail  to  illustrate  the  methods  of  calculation, 
consistent  with  the  requirement  of  reasonable  brevity,  will  be  presented. 

Consider  a  vector  rotated  through  an  angle  0  about  an  axis  n^  into 
the  vector  k,  lying  on  the  cone  shown  in  Figure  1 .  Simple  vector  analysis 
shows  that  the  vector  k^  A  (n^  A  )^)  is  normal  to  the  plane  tangent  to  the 
cone  and  containing  the  unit  vector  normal  to  this  plane  is  given  by 


S  -  £04^ 
(1  -  4-^)^ 


(C.1) 


where  k^ 


and  4^  -  ISo’— o  '  cosine  of  the  angle  between  k^,  n^. 

— o 


From  (18)  we  see  that 

k  -  k  |(  1  -  cos 6  )  ijjn  +  cos6  k^  +  sinOn^  A  k^j  ( C .  2  ) 

and  so 


§jL-]So  ■  ^4'(1  -COs6)(1 


(C.3) 


the  component  of  k  resolved  onto  the  tangent  plane  is  then  given  by 


liH  ■  is  ■  "  (4*^  +  C  -  4J^)cose)kQ  +  sin0n  A  Kq  (C.4) 


We  now  consider  the  inverse  of  this  problem;  given  a  vector  lying  in 
the  tangent  plane,  what  is  the  corresponding  g_j_so  that  q  =  +  qj_  lies  on  the 
cone,  qj^  is  resolved  into  components  within  the  tangent  plane  parallel  and 
perpendicular  to  k^  given  by 


Cl 


3H-iio 


(Ko  ^  Ho) 


(1  - 


These  equations  can  now  be  solved  for  0,  the  angle  which  rotates  qk-  into  g. 
From  (C.5)  we  see  that 


q(i()2  +  (1  -  4)2)cos0) 


-qsinS  (1  -  )  • 


Now  on  introducing  3/a  =  n  we  find  that  6  satisfies 


(C.6) 


+  COS0 


(C.7) 


where  ^  . 

To  solve  this  we  make  the  substitution  C  -  e^^  which  converts  (C.7)  into  the 
quadratic  equation 


c2(i  -  A)  .  2C  ^  ^ 


whose  solutions  are 


(1  -  S^)  ^ 

c2 


(1  -  3^)2  ( 

1  -  M 

'J  s4  i 

n^s^J 

We  see  that  as  long  as 


1  ^  (1  - 


(C.10) 


ICl -  1  and  6  is  real.  Geometrically  this  corresponds  to  the  condition  that 
IS  under  the  cone.  To  identify  the  appropriate  root  of  the  quadratic 
equation  we  let  g^  become  collinear  with  .  Then  3  -►  o,  q  o  and 


C2 


c  -  ±. 


TTns' 


for  C  -  1  as  6  -  0  we  chose  the  root 


+  i  1 


(1  -  5^)2 


(C.11) 


From  this  we  see  that 


0  =  Ti  +  4)^  -  4)2 


(C.12) 


where 


4>1  -  tan  ‘’I - ^  (s2  +  r)2)  -  i 

1(1  -  J 


(C.13) 


<J)2  -  tan" 


-1  /  1 


“  ti/2  <  4>i  /  4>2  ^  ti/2  . 


Having  found  6,  q  is  given  by 


{\\}^  +  ( 1  -  4^^)cos0) 


(C.14) 


so  that 


3.  -  (n„  - 

+  (1  -  ip^)cos0) 


(C.15) 


C3 


while  it  is  possible  to  express  cos  0  in  terms  of  4^,r{  the  resulting 
expressions  are  cumbersome  and  of  little  practical  use  and  consequently  are 
omitted. 

We  now  require  an  expression  analogous  to  (15)  for  close  to  . 
This  can  be  obtained  by  direct  expansion  of  the  results  we  have  just  obtained. 
In  outline,  we  set 


~  ^^1 

(3  =  6q2 

so  that 


and  find  to  lowest  order  in  6q^ ,  6q2 


It  is  reassuring  to  note  that  this  limiting  form  of  the  result  is  in  agreement 
with  that  obtained  in  the  previous  appendix  (eq  B.19) 
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The  vector  is  rotated  about  an  axis  n  into  k,  lying  on  the 
surface  of  a  cone. 


Figure  1 . 
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Figure  2(a) 
Simulation  Case  1 
Radar  and  target  geometry 

Target  coordinates: 
a*  (0,67,0),  b:  (67,0,0) 

c:  (67,57,0) 

Radar  parameters: 

Centre  frequency  fo  =  10  GHz, 
Target  rotation  rate  =  4®  S”^, 


The  radar  transmission  is  of  the  wi 
than  stepped  frequency* 


Figure  2(b) 
Simulation  Case  2 
Radar  and  target  geometry 


Target  coordinates: 


a :  (0,67,0),  b : 

(67,0,0) 

c:  (67,67,0),  d: 

(0,0, 133) 

Standard  deviation  of 

bandwidth: 

=  75  MHz 

Standard  deviation  of 
time 

integration 

=  125  mS. 

instantaneous  bandwidth  type,  rather 
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